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ABSTRACT 

We study the interaction between the atmospheres of Venus-like, non-magnetized exoplanets orbiting 
an M-dwarf star, and the stellar wind using a multi-species Magnetohydrodynaic (MHD) model. We 
focus our investigation on the effect of enhanced stellar wind and enhanced EUV flux as the planetary 
distance from the star decreases. Our simulations reveal different topologies of the planetary space 
environment for sub- and super-Alfvenic stellar wind conditions, which could lead to dynamic energy 
deposition in to the atmosphere during the transition along the planetary orbit. We find that the 
stellar wind penetration for non-magnetized planets is very deep, up to a few hundreds of kilometers. 

We estimate a lower limit for the atmospheric mass-loss rate and find that it is insignificant over the 
lifetime of the planet. However, we predict that when accounting for atmospheric ion acceleration, a 
significant amount of the planetary atmosphere could be eroded over the course of a billion years. 

Subject headings: planets and satellites: atmospheres - planets and satellites: magnetic fields - planets 
and satellites: terrestrial planets - magnetohydrodynamics (MHD) 


1. INTRODUCTION 

Nearly two decades after the first detections of exo¬ 
planets (Mayor & Queloz 1995; Schneider 1995; Mayor 
et al. 2003), and following the Kepler mission (Borucki 
et al. 2010), the focus of exoplanetary research has grown 
from individual planet detections to studies of the sta¬ 
tistical distributions of planets. In addition, the growth 
and quality of available data now enables us to study 
and characterize individual exoplanets in more quantita¬ 
tive detail. 

In the search for exoplanets, much interest is now di¬ 
rected to planets orbiting their host star in the Habit¬ 
able Zone (HZ). The commonly used definition of the 
HZ — the potential for a planet to have liquid water 
on its surface as a function of host stellar luminosity and 
planetary orbital separation (e.g., Kasting et al. 1993) — 
limits the physics involved in characterizing exoplanets 
to the radiation forcing applied at the top of a planetary 
atmosphere. Among other effects that can impact hab¬ 
itability are atmospheric loss due to intense EUV/X-ray 
radiation (focusing on hot jupiters, e.g., Lammer et al. 
2003; Baraffe et al. 2004; Yelle 2004; Baraffe et al. 2006; 
Garcia Munoz 2007; Murray-Clay et al. 2009; Adams 
2011; Trammell et al. 2011; Owen & Jackson 2012), cloud 
coverage and albedo, greenhouse gases, and atmospheric 
circulation (e.g., Tian et al. 2005; Cowan & Agol 2011; 
Heller et al. 2011; van Summeren et al. 2011; Wordsworth 
& Pierrehumbert 2013). 

Despite the ever more detailed characterization of the 
atmospheres of exoplanets, very few studies have inves¬ 
tigated their space environment and its possible physi- 
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cal impact. Of particular importance is the ability of a 
close-in planet to sustain an atmosphere at all. As an ex¬ 
ample from our own solar system, it is believed that the 
Martian atmosphere has been significantly eroded and 
evaporated over time (e.g., Lammer et al. 2013, with ref¬ 
erences therein). A number of studies have investigated 
the interaction between the atmospheres of close-in plan¬ 
ets and the surrounding plasma using hydrodynamic par¬ 
ticle models, where they used scaled stellar wind param¬ 
eters (Erkaev et al. 2005; Kislyakova et al. 2013, 2014), or 
scaled Coronal Mass Ejection (CME) parameters (Kho- 
dachenko et al. 2007; Lammer et al. 2007) as the input at 
the top of their models. All of these studies have found 
high escape rates of different forms from close-in planets, 
partially due to the intense EUV radiation. 

Sterenborg et al. (2011) employed detailed magneto¬ 
hydrodynamic (MHD) wind models tailored to the case 
of the young, active Sun to show that the magnetopause 
standoff distance of the Archean Earth could have been 
as small as 4.25 Earth radii (Re) in the face of a much 
more vigorous solar wind, compared with the present-day 
distance of 10.77?^; (see also Tarduno et al. 2010). Simi¬ 
larly compressed planetary magnetospheres might be ex¬ 
pected for exoplanets around other stars of high mag¬ 
netic activity (e.g., Khodachenko et al. 2007; Johansson 
et al. 2009, 2011a,b; Vidotto et al. 2013; See et al. 2014). 
More recently, (Cohen et al. 2014, Col4 hereafter), stud¬ 
ied the interaction of stellar winds with Earth-like, mag¬ 
netized exoplanets orbiting M-dwarf stars at habitable 
zone distances. They found that such close-in planets 
can transition between sub- and super-Alfvenic plasma 
environments along their orbits, and that this dynamic 
transition leads to a major heating at the top of the plan¬ 
etary atmosphere. 

The interaction of non-magnetized planets with a stel¬ 
lar wind is very different to such an interaction of mag¬ 
netized planets. In the former case, the stellar wind di¬ 
rectly interacts with the atmosphere and not with the 
planetary magnetic field and magnetosphere (Kivelson 
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& Russell 1995). In the solar system, such an interaction 
occurs at Mars, which has a very thin atmosphere and 
weak crustal field in its southern hemisphere (e.g., Ma 
et al. 2004a). Direct interaction between the solar wind 
and the atmosphere also occurs at Venus. In the Venu¬ 
sian case, the atmosphere is thick enough to balance the 
pressure from the solar wind (Hunten et al. 1983). 

Modeling the interaction of non-magnetized bodies 
with a flow of magnetized gas is very complicated, be¬ 
cause an accurate approach needs to account for the at¬ 
mospheric chemistry, ionization and radiation processes, 
and the atmospheric acceleration processes. All of these 
processes define the planetary mass source that resists 
the stellar wind plasma flow, and the induced magneto¬ 
sphere created as a result of this interaction. A num¬ 
ber of models have been developed to study this type of 
problem, including the interaction of the solar wind with 
Mars (e.g., Ma et al. 2004a; Ma & Nagy 2007, with ref¬ 
erences therein), and the interaction of Titan with the 
Kronian (Saturn) magnetosphere (e.g., Ma et al. 2004b, 
2006, with references therein). 

Ma et al. (2013) (Mal3 hereafter) have studied the in¬ 
teraction between the Venusian atmosphere and the so¬ 
lar wind. They focused on studying the shock location 
in front of the planet for solar minimum and solar max¬ 
imum conditions, and on the overall three-dimensional 
structure of the nearby space environment. Their model 
compared well with observations. 

In this paper, we use the model by MaI3 to study 
the interaction between Venus-like exoplanets around M 
dwarfs and the stellar wind of their host star. Since 
the parameters describing exoplanetary systems are cur¬ 
rently only loosely constrained, we leave many features 
of the model unchanged. Here, we focus on close-in plan¬ 
etary orbits and the competing effects of the enhanced 
upstream stellar wind and the increased EUV radiation. 
These are the two main changes we apply to the original 
modeling work of Mai 3. 

The structure of the paper is as follows. In the next 
section, we describe the numerical model and the over¬ 
laying assumptions we make. We describe the results in 
Section 3 and discuss their consequences in Section 4. 
We summarize our conclusions in Section 5. 


ing the local time stepping method (to accelerate conver¬ 
gence time Toth et al. 2012). In the case of sub-Alfvenic 
stellar wind conditions, we impose a float boundary con¬ 
dition for the pressure in order to mimic the two-way 
interaction between the incoming flow and the flow in 
the simulation domain. 

All the simulations presented here were performed on 
the Smithsonian Institute’s HYDRA cluster. 

2.2. MHD EQUATIONS 
The model solves the following set of equations. 

Multi-species continuity equations: 


^ + V • (ftu) = 

Si = niiNi (^lyph + E ksiUs^ , 


Li = m^Ui ( aR^iUe -t E kitUt ) , (1) 

t 


where the index i stands for the different ion species used 
for Venus - H+, 0+, 0^ and CO^ 

Here, pi are the mass densities of the different ions 
and the corresponding number densities are denoted by 
Ui, mi is the mass of the particular ion, and Ni is the 
number densities of the neutral sources, CO 2 and O. The 
ion source term, Si, depends on the photoionization rate, 
Vph,i^ the impact ionization rate, Vimp,i, and the charge 
exchange rate constant, kai (the index s sums over all 
the ions). The ion loss term, Li, depends on the charge 
exchange rate constant, ku (the index t sums over all 
the neutrals), and the recombination rate constant, a/jy. 
The rates and rate constants of the chemical reactions 
considered here are summarized in Table 1. 

Once the number densities of the different ions are 
determined, the total plasma density is defined as 
p = ^ Pi- This plasma density is then used in the rest 

i 

of the model equations. 


Momentum equation: 


2. NUMERICAL SIMULATIONS 
2.1. SIMULATIONS SETUP 

All the simulations presented here were done using the 
generic BATS-R-USMHD code (Powell et al. 1999; Toth 
et al. 2012), in a Cartesian geometry. The simulation do¬ 
main is defined as —40i?t, < X,Y, Z < 40Ry, with Venus 
radius, Rv = 6052 km. The simulation coordinate sys¬ 
tem is defined as follows. The X axis is defined as the line 
joining the centers of the planet and star, the Z axis is 
perpendicular to the ecliptic plain, and the Y axis com¬ 
pletes the right-hand system. We use initial mesh refine¬ 
ment to generate a grid structure with high resolution 
concentrated around the planet, with a grid cell size of 
60 km between the surface and an altitude of about 600 
km. This high resolution enables us to better resolve the 
induced ionosphere of the non-magnetized planet. The 
stellar wind upstream conditions are imposed as fixed 
boundary conditions at the box side facing the star and 
the simulation is run until a steady-state is achieved us¬ 


^-fV.(puu + pI+|^I-iBB) 

= pG - E P* E - E (2) 

it i 

Here, u is the (single fluid) plasma velocity vector, 
p is the pressure, B is the magnetic field vector, 
G is the gravitational acceleration, po is the mag¬ 
netic permeability in vacuum, I is a unit matrix, and 
t'ty = 4-10“^°{[0]-|-[C'02]} is the collision frequency 
between ion i and neutral t (Schunk & Nagy 1980). 

Induction equation: 

t+V-(uB-Bu)=Vx (^Vxh). (3) 

This equation describes the impact of the ambient 
plasma on the magnetic field, where the right-hand side 
term describes the effect of the induced magnetosphere, 
with conductivity cto- cto depends on the electron mass. 
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We, the electron number density, Ue, the electron charge, 
e, and the atmospheric collision frequencies between elec¬ 
trons and ions {vei), and between electrons and neutrals 
{yen)'- 


CTo = 


neC" 


me{Ver + J^en) ' 


( 4 ) 


The collision frequencies are taken from Schunk & 
Nagy (1980) and they both depend on the electron 
temperature, Tg. Vgi also depends on the electron 
density, rie, while t'en depends on the neutral densities, 
[O] and [CO 2 ]. 


Energy equation: 


— -I- V 

dt ^ ^ 




VxB 


= PU-G-J2J2 




, (5) 


where e = \pu^ + + 2 ^^^ is th® energy density, 

7 = 5/3 is the ratio of specific heats, nrit is the mass of 
neutral t, k is the Boltzmann constant, and T„ = IbOOiti 
is the temperature of the newly created ions. In this 
model we assume that Tg and the ion temperature, T^, 
are the same and that they are both equal to half of the 
plasma temperature, Tp = kp/nipp. 


2.3. STELLAR WIND INPUT 

For the upstream stellar wind conditions of planets or¬ 
biting an M-dwarf star, we follow the work presented in 
Col4. The stellar wind conditions of EV Lac, an M3.5 V 
star, were obtained using an MHD model for the stel¬ 
lar corona and wind. This model was driven by available 
surface magnetic held data of the star (Morin et al. 2008, 
see Col4 for more details on the coronal model). Once 
the stellar wind solution is obtained, the particular con¬ 
ditions are extracted at the orbits of three potentially 
habitable planets observed by Kepler - Kepler Object 
of Interest (KOI) 2626.01, 1422.02, and 854.01 (Dress¬ 
ing & Charbonneau 2013). These planets are located 
at distances of 0.06, 0.085, and 0.15 AU, respectively. 
For each of these planets, we study two cases in which 
the planets reside in sub- or super-Alfvenic stellar wind 
conditions. For reference, we simulate Venus itself at 
an orbit of 0.72 AU and a ’’Venus at the orbit of Mer¬ 
cury” case at 0.39 AU. For these two cases, the stellar 
wind is super-Alfvenic, but we study a case of fast and a 
case of slow stellar wind. We also assume here that the 
interplanetary magnetic held (IMF) has a pure Bz com¬ 
ponent. This is due to the fact that the planet does not 
have an internal held, and that the IMF direction does 
not seem to affect the O’*" escape from Venus (Masunaga 
et al. 2013). The upstream stellar wind conditions used 
for all cases are summarized in Tables 2 and 3. 


2.4. ENHANCED STELLAR EUV ELUX 

In MaI3, the changes in EUV radiation over the solar 
cycle were considered with two cases for solar minimum 


and solar maximum. Their model also included the ra¬ 
diation change with Solar Zenith Angle (SZA). 

Planets that reside in close-in orbits are expected to 
receive enhanced amounts of X-ray and EUV radiation 
compared to the Earth, Venus, and even Mercury. In or¬ 
der to consider the effect of this enhanced EUV radiation, 
we scale the photoionization rates of the solar maximum 
case from MaI3 with the distance square. The scaling pa¬ 
rameter of Venus’s EUV flux, Fc^, for the different cases 
is shown in the last column of Tables 2 and 3. 

We stress that here we expect the enhanced EUV ra¬ 
diation to simply increase the creation of ions via pho¬ 
toionization processes. In our simplified model for the at¬ 
mosphere, we do not account for the overall atmospheric 
heating due to the close proximity to the star (since we 
are considering the habitable zones of M dwarfs, the op¬ 
tical and IR heating will be similar to the venus-earth- 
mars case), nor to the possible change in the electron 
temperature as a result of the creation of super-thermal 
electrons (i.e., we do not scale Te, Ti, and Tn with the 
EUV flux). These effects can potentially lead to charge 
separation as the result of the change in ion and electron 
scale heights, and to additional acceleration and escape 
of atmospheric ions (i.e., the polar wind, Cohen & Glocer 
2012). We leave the implementations of these processes 
in the model for future work. 

2.5. ATMOSPHERIC CHEMISTRY, THE 
NEUTRAL ATMOSPHERE, AND INNER 
BOUNDARY CONDITIONS 

The model assumes a neutral planetary atmosphere 
source composed of Carbon dioxide, which dominates be¬ 
low 160 km, and Oxygen, which dominates above that 
height (at Venus). The atmospheric neutral densities as¬ 
sumed here are (the solar maximum case in Mai 3, Fox 
& Sung 2001): 

[CO 2 ] = I.O • lO^^ • cto-3 

[O] = 2.0 • lO^i • cm-3, ( 6 ) 

where the altitude is in km, and zq = 100 km is the inner 
boundary of the model. 

As seen from Eq. I, we assume that the atmospheric 
chemistry is driven by these two neutral species, and it 
is governed by the photoionization, impact ionization, 
charge-exchange, and recombination processes described 
in Table I. In future work, we plan to implement other 
species in order to study exoplanetary atmospheric envi¬ 
ronments. 

The densities for the ion species at the inner boundary 
are specified by the photochemical equilibrium. The in¬ 
ner boundary is set to absorb the velocity and magnetic 
field vectors. This guarantees that the magnetic field 
diffuses through the boundary with the help of the low 
conductivity (as a result of the strong ion-neutral cou¬ 
pling at these low altitudes). The plasma temperature is 
assumed to be twice the (Venusian) neutral temperature 
of 200 K. Here we also neglect a potential increase of the 
neutral temperature due to close proximity to the star 
(assuming that such an increase is possible considering 
the actual stellar radiation of M-dwarfs). 
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3. RESULTS 

Figure 1 shows the model solutions in the Y=0 and 
Z=0 planes for all cases. The star is on the right, and 
the colored contours indicate the ratio of ionized Oxygen 
(0++0^) to H+ densities. The planet is represented by 
a sphere and selected magnetic field lines and also shown. 

It can be seen that in all the cases, parts of the plane¬ 
tary environment are dominated by the oxygen ions over 
the hydrogen ions by factors of a few tenths (and up to 
a thousand near the inner boundary). The exception is 
the case for 0.06 AU with a super-Alfvenic wind, which 
remains dominated by the stellar wind hydrogen popula¬ 
tion. In all other cases, the lower parts of the atmosphere 
are completely occupied by oxygen, as are large regions 
in the planetary wake (these are occasionally referred to 
as the ’’induced magetosphere” and ’’induced magneto¬ 
tail” ). 

A notable feature is that the oxygen-dominated regions 
are largely confined to the equatorial plain in most of 
the super-Alfvenic cases, while these regions extend to 
lobes above and below the planet in the sub-Alfvenic 
case. The super-Alfvenic Mercury case also shows similar 
structure due to the low Alfvehic Mach number (close to 
2) in this case. The structure that eappears in the sub- 
Alfvenic cases is typical for a body that moves in ambient 
plasma with sub-Alfvenic speed, where the ambient field 
is perturbed by the body via a unique current system 
known as ’’Alfven Wings” (Neubauer 1980, 1998). In 
this topology, which has been studied in the context of 
the Jovian moons of lo and Ganymede (Combi et al. 
1998; Linker et al. 1998; Jacobsen et al. 2007; Ip & Kopp 
2002; Kopp & Ip 2002; Jia et al. 2008), two standing 
“lobes” are created above and below the body. In our 
simulations, these lobes are filled with the oxygen ions 
that originate from the planetary atmosphere. 

Figure 2 shows the velocity distribution for all the cases 
on the Y=0 (meridional plain) and the Z=0 (equato¬ 
rial plain), together with the contour line at which the 
Alfvenic Mach number equals one. One can see that 
there is no bow shock in front of the planet for the sub- 
Alfvenic cases, and that the stellar wind is directly slowed 
down to speeds of less than 100 km/s at a altitude of 
about 600-800 km at the sub-stellar point (it is even lower 
in altitude at the planetary flanks). In the super-Alfvenic 
cases, a shock is formed at altitudes of at least twice this 
height, but the stellar wind incoming speeds are rather 
high (about 100 km/s) at much lower altitudes. 

Another representation of the Alfven wing topology 
is the vertically extended planetary wake with very low 
speeds seen in Figure 2 Y=0 plains for the sub-Alfvenic 
cases and the super-Alfvenic cases with low Alfvenic 
Mach number (0.06AU and Mercury). The super- 
Alfvenic cases with strong shock (i.e., high Alfvenic Mach 
number) show a wake that is much more conhned towards 
the equator. Overall, all cases show regions around the 
planet with a very low flow speed, which indicates a stel¬ 
lar wind free zone. 

The velocity pattern in the sub-Alfvenic cases shows 
that the flow is deflected to very large angles as the re¬ 
sult of the standing vertical lobes. In the Super-Alfvenic 
cases, the flow simply goes around the planet, generating 
a slightly turbulent wake behind the planet in the slow 
speed tail region. 


Figure 3 shows the total ion density, and the densities 
of the different ion species, as a function of altitude for 
all the cases. The densities are extracted from the plane¬ 
tary surface along the sub-stellar line (X axis). The total 
ion density is dominated by H+ down to an altitude of 
1000-1500 km, while below this height the atmospheric 
ions begin to dominate the ion density and peak around 
300 km. The enhanced creation of atmospheric ions as 
a result of the increased EUV flux is clearly reflected in 
an increase in ion densities as a function of decreasing 
orbital distance. The plots also show that the altitudes 
at which the densities of the atmospheric ions are com¬ 
parable with the density of stellar wind H+ are located 
lower in the atmosphere for the super-Alfvenic conditions 
compared to the sub-Alfvenic conditions (for the three 
close-in cases). This is due to the fact that the stellar 
wind density is much higher for the super-Alfvenic cases. 

4. DISCUSSION 

The results of our simulations show that for close-in 
Venus-like exoplanets, a very strong, direct interaction 
between the stellar wind and the upper atmosphere oc¬ 
curs at an altitude of no more than about one planetary 
radius (a few thousand kilometers). This is in contrast to 
the few planetary radii for the case of magnetized planets 
presented in Gold. 

4.1. TOPOLOGY OF THE NEARBY SPACE 
ENVIRONMENT 

The results show a quite distinct topology of the 
plasma environment near the planet for the sub- and 
super-Alfvenic stellar wind upstream conditions. This 
suggests that the space environment of close-in Venus- 
like planet can be very dynamic as these planets pass be¬ 
tween sub- to super-Alfvenic plasma environments along 
their orbit (as demonstrated dynamically in Col4 but ex¬ 
cluded here). Therefore, through planetary orbital mo¬ 
tion alone (regardless of any stellar activity), the stellar 
wind can affect the planetary atmosphere in terms of dif¬ 
ferent forms of energy deposition, such as Ohmic dissipa¬ 
tion, Joule heating, and the generation of gravity waves. 
All these are expected due to the dynamic change of the 
induced ionosphere-magnetosphere. 

4.2. STELLAR WIND PENETRATION 

Figure 4 shows the plasma radial speed as a function 
of altitude up to about 3 planetary radii for the magne¬ 
tized planets from Col4 and the non-magnetized planets 
studied here. The plot shows the cases corresponding to 
an orbital distance of 0.06 AU. A negative sign of this 
speed component represents an inflow, while a positive 
sign represents an outflow. The plot clearly shows that 
for the magnetized planet cases, the inflow is already very 
weak far from the planet and there is even weak outflow 
in the sub-Alfvenic case. In the case of non-magnetized 
planets, the inflow slowly decreases for the super-Alfvenic 
case (this is the post-shock flow), and sharply decreases 
for the sub-Alfvenic case. 

Figure 5 shows the radial speed and the magnitude 
of the magnetic field as a function of altitude for all 
the non-magnetized Venus-like cases. The values are ex¬ 
tracted from the planetary surface along the sub-stellar 
line (similar to Figure 3). The scale for the radial speed is 
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designed to demonstrate the point at which the incoming 
wind gets to a complete stop (inflow speed of less than 
1 km/s). In both the sub- and super-Alfvenic cases, line 
1 marks the point where the wind stops in the cases of 
Venus, and line 2 marks the point where the wind stops 
in the case of the planet located at 0.06 AU. In both 
cases, the altitude difference is about 200 km and the 
wind penetrates up to an altitude of less than 600 km 
above the planetary surface. 

The magnetic held slowly increases as the result of the 
plasma compression in front of the planet, until a max¬ 
imum is attained, and then starts to decrease to zero 
at an altitude of around 180 km for the case of dis¬ 
tance of 0.06 AU (marked by line 3). This point, which 
marks the decline of the magnetic held and represents 
the magnetic field-free planetary atmosphere is the same 
for all the cases of super-Alfvenic wind conditions (with 
the exception of the Mercury case which is almost sub- 
Alfvenic). In the sub-Alfvenic cases, this point moves 
outward as the planet gets further from the star, up to 
about 350 km for the Mercury and Venus cases. This dif¬ 
ference is probably due to the fact that the stellar wind 
in the sub-Alfvenic cases is much less dense than the 
super-Alfvenic one. Therefore, the compression region of 
the incoming stellar wind held, as well as the ability of 
the field to diffuse towards the atmosphere (as the result 
of increased ion/electron density and their effect on the 
conductance and collision frequencies), are more sensi¬ 
tive to the Oxygen outflow in the sub-Alfvenic case than 
the super-Alfvenic case. 

To demonstrate the strong effect of the increased EUV 
radiation, Figure 6 shows the difference in the radial 
speed as a function of altitude for the case located at 
0.06 AU with and without scaling of the EUV radiation 
with distance. The plot clearly shows that the enhanced 
radiation increases the production of atmospheric ion 
flux, which stops the stellar wind about 200 km higher 
than the case without increase EUV flux. This increased 
ion flux also results in an enhanced compression of the 
magnetic held, even though it does not seem to affect the 
point at which the magnetic field starts to decrease. 

4.3. ATMOSPHERIC ION ESCAPE 

The model used here separates the densities of the dif¬ 
ferent ions in the continuity equation but it does not dis¬ 
tinguish between the different ion velocities and uses a 
single-fluid plasma velocity in the momentum, induction, 
and energy equation. Since we do not have information 
about the particular ion velocities, it is hard to estimate 
the escape rate of the different ions, in particular the 0+ 
escape rate. In addition and as mentioned in Section 2.4, 
the particular model used here does not account for any 
acceleration of the atmospheric ions via hydrodynamic 
acceleration (e.g., Garcia Munoz 2007; Tian et al. 2005), 
and the ambipolar electric field (Axford 1968; Banks & 
Holzer 1968; Cohen & Glocer 2012), which is the main 
mechanism to explain the large amount of O'*" observed 
in the Earth’s magnetosphere. Tian et al. (2009); Tian 
(2009) even found that the intense EUV radiation in¬ 
creases the amount of ionized CO 2 and decrease the 
amount of neutral CO 2 . As a result, the atmospheric In¬ 
frared (IR) cooling is reduced as well and the atmosphere 
becomes over-inflated. However, this is a hydrodynamic, 
one-dimensional model that does not account for day- 


night energy transfer, electromagnetic effects, and it does 
not include the dynamic pressure of the incoming stellar 
wind at the top boundary. 

Nevertheless, we can attempt to estimate a lower limit 
for the 0+ escape rate by calculating the the integral of 
the mass flux through a sphere at a certain height, and 
include only the density points at which the radial veloc¬ 
ity has both positive sign (i.e., outflow), and magnitude 
larger than Venus’s escape velocity of 10.36 km/s. Fig¬ 
ure 7 shows the loss rates of O'*' for the different cases, 
integrated on spheres of r = 300 km, which is the alti¬ 
tude of the ion density peak, and r = 600 km, which is 
twice that height. The values are not always the same 
for the two heights as our assumption here does not con¬ 
serves mass. In addition, some of the integrations result 
in zero as they do not include any outflow points with 
speeds above the escape velocity. 

The loss rates for the Venus case are of the order of 
10 ^^ and they are comparable with O’*" loss rates 
observed by Pioneer Venus Orbiter and Venus Express 
(e.g., hammer et al. 2006; Luhmann et al. 2008; Du¬ 
binin et al. 2011; Fedorov et al. 2011; Lundin 2011). 
The escape rate is enhanced by 3-4 orders of magni¬ 
tude for the close-in planets as the result of the enhanced 
EUV photoionization. The mass loss rates translate to 
10“^® — 10“^^ Mya yr~^, with Mya = 5 • 10^^ g be¬ 
ing the mass of the Venusian atmosphere (Basilevsky 
& Head 2003). Interestingly, the 0+ loss rates of the 
close-in planets are similar to those obtained for Venus 
4 Gyr ago by Kulikov et al. (2006), who calculated the 
Venusian loss rate of 0+ for different ages of the solar 
system (taking into account the changes in the stellar 
radiation environment). The loss rates obtained by our 
simulations are much lower than the loss rates calculated 
for atmospheric hydrodynamic escape from hot jupiters 
(10^ — 10^^ g e.g., hammer et al. 2003; Erkaev et al. 
2005; Penz et al. 2008; Murray-Clay et al. 2009; Owen 
& Jackson 2012). The reason is that hot jupiters lose 
hydrogen atoms via hydrodynamic blow off. In the case 
of stellar wind ion pick up (e.g., Kislyakova et al. 2013, 
2014, for the Kepler-11 planets), hydrogen-dominated at¬ 
mospheres are much more expanded compared to CO 2 
atmospheres, so that the pressure balance between the 
stellar wind plasma and the hydrodynamically expand¬ 
ing upper atmosphere is much further up compared to 
the cases presented here. 

In the case of a planet located at 0.06 AU, the loss 
rate is much higher for the sub-Alfvenic case comparing 
to the super-Alfvenic case. This is because the strong, 
dense stellar wind in the super-Alfvenic case suppresses 
the atmospheric Oxygen outflow generated by the en¬ 
hanced EUV radiation (see top-right panel of Figure 1). 
The outflow is much more significant for the sub-Alfvenic 
case, at which the stellar wind is less dense. For the cases 
of planets located at 0.085 AU and 0.15 AU, the loss 
rate is higher for the super-Alfvenic cases comparing to 
the sub-Alfvenic cases due to the slowdown of the stel¬ 
lar wind at the shock in front of the planet. This leaves 
more room for the atmospheric ions to outflow against 
the incoming wind. 

At face value, the loss rates estimated by our simula¬ 
tions mean that only a very small fraction of a Venusian- 
like atmosphere will be lost over the course of a billion 
years. Therefore, our simulations suggest that there is a 
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good chance that the atmospheres of Venus-like exoplan¬ 
ets orbiting M-dwarf stars are sustainable. Nevertheless, 
the actual mass-loss rate may be significantly higher than 
these lower limits. The increase in ion scale-height due to 
increased EUV photoionization may boost the ion den¬ 
sity by an order of magnitude at high altitude. In addi¬ 
tion, acceleration of ions at the upper atmosphere as a 
result of the polar wind can increase the ion velocity by 
another factor of 10 (Glocer et al. 2009; Cohen & Glocer 
2012). This would mean instead that the planet may lose 
most of its atmosphere over the course of a billion years. 
A similar conclusion was drawn by Lammer et al. (2007), 
who found that close-in terrestrial exoplanets may loose 
hundreds or even a thousand bars of their atmosphere 
due to the intense EUV radiation and CME atmospheric 
erosion, even in the case where the planet has a strong 
internal magnetic field. 

More rigorous estimates of the atmospheric loss would 
require coupling of the model for the planetary space en¬ 
vironment presented here with a model accounting for 
the structure of the upper atmosphere, that includes the 
acceleration of the atmospheric ions, the impact of the 
increased EUV on the atmospheric temperatures and the 
atmospheric conductance, and a self-consistent photoion¬ 
ization model. Such implementations are needed in order 
to provide a more accurate assessment of the ability of 
close-in exoplanets to sustain their atmospheres. 

5. SUMMARY AND CONGLUSIONS 

We study the interaction between the atmospheres of 
Venus-like, non-magnetized exoplanets orbiting M-dwarf 
stars and the stellar wind using a multi-species MHD 
model. We focus our investigation on the effect of en¬ 
hanced stellar wind and enhanced EUV flux as the plan¬ 
etary distance from the star decreases from Venus’s orbit 
at 0.72 AU to close-in orbits located at 0.15, 0.085, and 
0.06 AU. 

For the close-in orbits, we study the interaction be¬ 
tween sub- and super-Alfvenic stellar wind conditions. 
We find that the two cases result in a very different 
configuration of the immediate planetary space environ¬ 
ment, with the sub-Alfvenic case resulting in vertically 
extended wake, while the super-Alfvenic case produces 
a wake which is confined to the equatorial region. We 
expect that such dynamic change along the orbit of few 
days would deposit energy into the planetary atmosphere 
in the form of heating and the generation of gravity 
waves. 

The stellar wind penetration for the non-magnetized 
planets studied here is much deeper than the magnetized 
planets studied in Gold. It reaches altitudes of several 
hundreds of kilometers, in contrast to thousands in mag¬ 
netized planets. As the planets reside closer to the star, 
the increase EUV photoionization leads to an enhanced 
creation of ions that push the stellar wind stagnation 
point up by about 200 km compared to a low photoion¬ 
ization rate case. It also increases the compression of the 
stellar wind magnetic field in front of the planet. 

We estimate a lower limit for the O'*' escape rate and 
find that the planetary atmosphere may be sustainable 
over the lifetime of the planet. However, the atmospheric 
escape rate could be much higher when additional accel¬ 
eration mechanisms for the atmospheric ions are consid¬ 
ered. We plan to implement these mechanisms in future 


work in order to better estimate the ability of Venus-like 
exoplanets to sustain their atmospheres. 
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TABLE 1 

Chemical Reactions Used in the Model 


Reaction 

Rate CeofRcient 

Reference 

CO 2 + hiz —i- CO 2 + e 

3.24 X 10“® s-i 

Schunk Sz Nagy (2009) 

0 + hiz —i- 0 + + e 

1.21 X 10“® s-i 

Schunk & Nagy (2009) 

cot + 0 -t ot + CO 
Cof + 0 -t 0 + + CO 2 

1.64 X IQ-lO cm-® 

Schunk & Nagy (2009) 

9.6 X 10“®® cm”® s”® 

Schunk & Nagy (2009) 

0 + + CO 2 -t ot + CO 

1.1 X 10“® cm-® S-® for Ti < 800 K; 

1.1 X 10-®(800/ri)® ®® cm-® S-® for Ti > 800 K 

Fox & Sung (2001) 

H+ + 0 -t 0 + +H 

3.75 X 10-1® cm-® 5 -® 

Schunk & Nagy (2009) 

Ot + e 0 + 0 

1.95 X 10-^(300/re)0'^ cm-® s-® for Te < 1200 K-, 
7.38 X 10-®(1200/re)0 ®® cm-® 5 -® for Te > 1200 K 

Schunk & Nagy (2009) 

cot + e -t CO + 0 

3.5 X 10-’’(300/Te)® ® cm-® 5 -® 

Fox & Sung (2001) 


^ Electron impact ionization is neglected in the calculation, density is from the stellar wind, The neutral 

hydrogen is neglected. 


TABLE 2 

Sub-Alfvenic Stellar Wind Parameters 


r [AU] n 

[cm 

■®] T [K] 

[km/s] 

Uy [km/s] 

Uz [km/s] 

Bx [nT] 

By [nT] 

Bx [nT] 

Ma 

Feuv [r$] 

0.06 (KOI 2626.01) 

450 

1, 000, 000 

-600 

0 

0 

0 

0 

2,000 

0.29 

144 

0.085 (KOI 1414.02) 

160 

750, 000 

-700 

0 

0 

0 

0 

800 

0.5 

72 

0.15 (KOI 854.01) 

45 

500, 000 

-700 

0 

0 

0 

0 

250 

0.86 

23 

*0.39 (Mercury) 

15 

300, 000 

-700 

0 

0 

0 

0 

50 

2.5 

3.4 

*0.72 (Venus) 

10 

250,000 

-700 

0 

0 

0 

0 

10 

10.1 

1 

* 

The cases labeled ’’Venus” 

and 

’’Mercury” are super- 

-Alfvenic. 









TABLE 3 

Super-Alfvenic Stellar Wind Parameters 

r [AU] 

n [cm 

T[K] 

Ux [km/s] 

Uy [km/s] Uz 

[km/s] 

Bx [nT] 

By [nT] 

Bz [nT] 

Ma 

Feuv [-^9] 

0.06 (KOI 2626.01) 

35,000 

800, 000 

-200 

0 

0 

0 

0 

500 

3.4 

144 

0.085 (KOI 1414.02) 

12,500 

450, 000 

-250 

0 

0 

0 

0 

250 

5.1 

72 

0.15 (KOI 854.01) 

3,000 

200, 000 

-300 

0 

0 

0 

0 

100 

7.5 

23 

0.39 (Mercury) 

50 

300, 000 

-450 

0 

0 

0 

0 

50 

2.9 

3.4 

0.72 (Venus) 

10 

250,000 

-450 

0 

0 

0 

0 

10 

6.5 

1 
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Sub-Alfvenic Super-Alfvenic 



Fig. 1.— The solutions for all the cases shown on the Y=0 and Z=0 plains. The view is from the side where the stellar wind is coming 
from the right (the direction to the star). The left column shows the sub-Alfvenic cases and the right column shows the super-Alfvenic 
cases; the cases change from the closest planet (top) to the farther (Venus, bottom). Color contours are of the ratio of the Oxygen ion 
density (0"*"-|-Oj) to the H"*" density. Selected magnetic field lines are shown as white lines. 
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Fig. 2. — The plasma flow pattern for all the cases. Each pair of panel shows the meridional (X-Z) plain and the equatorial (X-Y) plain. 
The star is on the right. Color contours show the plasma velocity magnitude and the black lines show streamlines. The solid white line 
represents the contour for Ma = 1. 
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Sub- 

Alfvenic 




Super- 

Alfvenic 



Fig. 3. — Atmospheric density structure for the sub-Alfvenic (top) and the super-Alfvenic (bottom) cases extracted from the planetary 
surface along the sub-stellar line (X axis). Each panel shows the total density (solid), density (dashed), and ion density (dotted-dashed) 
as a function of altitude. The columns are for the ions of O'*" (left), Oj (middle), and COj (right). 
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R[km] 


Fig. 4. — Radial velocity as a function of altitude for the nearest non-magnetized cases calculated here (solid) and for the nearest 
magnetized cases calculated in Col4 (dashed). All cases located at 0.06 AU. Black lines represent sub-Alfvenic cases and red lines represent 
super-Alfvenic cases. 


Sub-Alfvenic Super-Alfvenic 



0.06 AU 
0.085 AU 

- 0.15AU 

Mercury 

Venus 

- - - 0.06 AU 

- - - 0.085 AU 

- - - 0.15AU 

- — — Mercury 

- — — Venus 


Fig. 5. — Magnetic field magnitude (solid lines) and radial velocity (dashed lines) as a function of altitude for all cases are extracted 
along the sub-stellar line (X axis). 
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Fig. 6 . — Magnetic field magnitude (left) and radial velocity (right) as a function of altitude for the nearest cases at 0.06 AU are extracted 
along the sub-stellar line (X axis). Black lines represent sub-Alfvenic cases and red lines represent super-Alfvenic cases. Solid lines represent 
the cases with EUV fiux scaled with distance while dashed lines represent cases with unsealed EUV fiux (Venus EUV fiux). 
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Fig. 7.— 0+ escape rate for the different cases as a function of distance from the star. The integration at 300 km did not yield any flux 
for the 0.06 AU, 0.15 AU, and Mercury cases. 






























